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^ I We have made the first measurement of the D*"*" width using 9/fb of e^e~ data collected 

■ near the T(4S) resonance by the CLEO II. V detector. Our method uses advanced tracking 

CO ' techniques and a reconstruction method that takes advantage of the small vertical size of the 

, CESR beam spot to measure the energy release distribution from the D*'^ — > D'^tt'^ decay. 

^ . Our preliminary result is r(D*+) = 96 ± 4 (Statistical) ± 22 (Systematic) keV. 

^ ■ 

o ; 

I A measurement of T{D*~^) opens an important window on the non-perturbative strong 

Q ■ physics involving heavy quarks. The basic framework of the theory is well understood, however, 

^ • there is still much speculation - predictions for the width range from 15keV to 150kevS. We 

QJ \ know the D*~^ width is dominated by strong decays, since the measured magnetic-dipole tran- 

' sition rate is small, Br{D*^ — > D'^j) = 1.68 it 0.45%, i and can be neglected to the first order. 

The level splitting in the B sector is not large enough to allow real strong transitions. Therefore, 
^ , a measurement of the £>*^ width gives unique information about the strong coupling constant 

^ \ in heavy-light systemsH, gD*D-K or g. 

Prior to this measurement, the D*^ width was limited to be less than 131 keV at the 90% 
confidence level by the ACCMOR collaborationB. The limit was based on 110 signal events re- 
constructed in two decay channels with a background of 15% of the signal. This contribution 
describes a measurement of the D*~^ width with the CLEO II. V detector where the signal, in ex- 
cess of 11,500 events, is reconstructed through a single, well-measured sequence, D*^ —>■ tt^q^-D'', 
— > K~TT^. Consideration of charge conjugated modes are implied throughout this paper. 
The level of background under the signal is less than 3% in our loosest selection. 

The CLEO detector has been described in detail elsewhere. All of the data used in this 
analysis are taken with the detector in its II. V configurationH. The data were taken in symmetric 
e^e~ collisions at a center of mass energy around 10 GeV with an integrated luminosity of 9.0/fb 
provided by the Cornell Electron-positron Storage Ring (CESR). The nominal sample follows 
the selection of D*^ — > tt^q^-D'^ "''"aiow-^"'""^ candidates used in our — mixing analysisB 
Our reconstruction method takes advantage of the small CESR beam spot and the kinematics 
and topology of the D*^ — > t^^iow-^^ ~^ '^tiow^^'^^ decay chain. The K~ and 7r+ are required 
to form a common vertex. The resultant candidate momentum vector is then projected 
back to the CESR luminous region to determine the production point. The CESR luminous 
region has a Gaussian width ~ 10 /im vertically and ~ 300 fim horizontally. This procedure 
determines an accurate production point. Then the vrj^^^ track is re-fit constraining its 



trajectory to intersect the production point. This improves the resolution on the energy 
release, Q, by more than 30% over simply forming the appropriate invariant masses of the 
tracks. The improvement to resolution is essential to our measurement of the width of the D*^. 
Our resolution is shown in Figure |l| and is typically 150 keV. The good agreement reflects that 
the kinematics and sources of errors on the tracks, such as number of hits and the effects of 
multiple scattering in detector material, are well modeled. 

The challenge of measuring of the width is understanding the tracking system response 
function since the experimental resolution exceeds the width we are trying to measure. We 
depend on exhaustive comparisons between a GEANT Q based detector simulation and our 
data. We addressed the problem by selecting samples of candidate D*^ decays using three 
strategies. First we produced the largest sample from data and simulation by imposing only 
basic tracking consistency requirements. We call this the nominal sample. Second we refine 
the nominal sample selecting candidates with the best measured tracks by making very tight 
cuts on tracking parameters. There is special emphasis on choosing those tracks that are well 
measured in our silicon vertex detector. This reduces our nominal sample by a factor of thirty 
and, according to our simulation, has negligible contribution from tracking mishaps. We call 
this the tracking selected sample. A third alternative is to select our data on specific kinematic 
properties of the D*~^ decay that minimize the dependence of the width of the D*^ on detector 
mismeasurements. The nominal sample size is reduced by a factor of three and a half and, again 
according to our simulation, the effect of tracking problems is reduced to negligible levels. We 
call this the kinematic selected sample. Table || summarizes the statistics in our three samples. 

In all three samples the width is extracted with an unbinned maximum likelihood fit to the 
energy release, Q = m{K~ 7:~^7r^^^) — m{K~7T~^) — m^+, distribution and compared with the 
simulation's generated value to determine a bias which is then applied to the data. The three 
different samples yield consistent values for the width of the D*~^ giving us confidence that our 
simulation accurately models our data. 

We assume that the intrinsic width of the is neg ligible, r(L>0) < F (£>*+), implying that 
the observed width of Q distribution is dominated by the shape given by the D*^ intrinsic 
width and the tracking system response function. Thus we consider the pairs of Qi and agi for 
D*~^ — > T^^io^D^ — > '''"Jow-^"'""^ where aqi is given for each i-th. candidate by propagating the 
tracking errors in the kinematic fit of the charged tracks!. We perform an unbinned maximum 
likelihood fit to the Q distribution, minimizing the likelihood function 

N 

L = 2{Ns + Nb)-2J2 ^og[Ns ■ SiQi,aQi; Fq, Qo, fmis, c^™,) + Nb ■ B{Q,- 61,2,3)], (1) 

i=l 

where S and B are respectively the signal and the background shapes, A'^, and number of 
signal and background events. 

The shape of the underlying signal is assumed to be given by a P-wave Breit-Wigner, with 
central value of Q, Qq. We considered a relativistic and a non-relativistic Breit-Wigners as a 
model of the underlying signal shape, and found negligible difference in the fit parameters. The 
width of the Breit-Wigner depends on Q and is given by 

where Fq is equivalent to F(D*+), m and P are the measured candidate D*^ mass and tt^^^ or 
momentum in the D*~^ rest frame and Pq and mo are the values computed using Qq. The 
effect of the mass term is negligible at our energy. The partial width and the total width differ 
negligibly in their dependence on Q for Q > 1 MeV . We use Eqn Q suitably normalized to 
describe of Q effects on width. 



(2) 



Table 1: Summary of our data sample, simulation biases, and fit results. 



Parameter 


Sample 


Nominal 


Tracking 


Kinematic 


Candidates 
Background Fraction (%) 

Tfit — Tgenerated (keV) 

Fit Fo (keV) 
D*+ Width (keV) 


11496 
2.51 ±0.27 
2.7 ±2.1 
98.9 ± 4.0 
96.2 ±4.0 


368 
4.1 ± 1.9 
1.7 ±6.4 
106.0 ± 19.6 
104 ± 20 


3284 
4.05 ± 0.49 
4.3 ±3.1 
108.1 ± 5.9 
103.8 ±5.9 



Table 2: Results of the fits described in the text. The uncertainties are statistical. 



Parameter 


Sample 


Nominal 


Tracking 


Kinematic 


Fo (keV) 


98.9 ± 4.0 


106.0 ± 19.6 


108.1 ± 5.9 


Qo (keV) 


5853 ± 2 


5854 ± 10 


5850 ± 4 




11207 ± 109 


353 ± 20 


3151 ± 57 


fmis (%) 


5.3 ±0.5 


NA 


NA 


amis (keV) 


508 ± 39 


NA 


NA 


Nb 


289 ± 31 


15 ±7 


133 ± 16 



For each candidate the signal shape, S, is a convolution of the Breit-Wigner function with 
a resolution Gaussian with width, crgj, determined by the tracking errors, as a model of our 
finite resolution. Figure |l| shows the distribution of aq for the data and the simulation. We 
allow a small fraction of the signal, fmis , to be parameterized by a single Gaussian with effective 
resolution, amis, different from measured, aqi. This shape is included in the fit to model the 
tracking mishaps, mis-assigned hits and hard multiple scatters, which our simulation predicts 
to be at the 5% level in the nominal sample and negligible in both the tracking and kinematic 
selected samples. For the purpose of systematic study the aq has a scale factor, k, which is 
fixed to one in our nominal fits. 

The fit also includes a background contribution with fixed shape, B, presented by polynomial 
function with three or more parameters, ^'1,2,3- This shape is taken from fits to the background 
prediction of our simulation. The level of the background is allowed to float in our standard fit. 

The fitter has been extensively tested both numerically and with input from our full sim- 
ulation. We find that the fitter performs reliably giving normal distributions for the floating 
parameters and their uncertainties. It also reproduces the input T(D*^) from to 130 keV with 
offset consistent with zero, as shown in Table |l[ 

We note that if all the parameters are allowed to vary simultaneously there is strong corre- 
lation among the intrinsic width Fq, the fraction of mismeasured events fmis, and the aq scale 
factor k, as one would expect. Thus our nominal fit holds k fixed to one, but in our systematic 
studies we either fix one of the three or provide a constraint with a contribution to the likelihood 
if the parameter varies from its nominal value. 

Figures ^, ^, and |^ respectively display the fit to the nominal, tracking, and kinematic 
selected data samples. The results of the fits are summarized in Table |2|. Correlations among 
the floating parameters of the fit are negligible. 

Figure |5| displays the likelihood as a function of the width of the D*^ for the fits to the three 
data samples. 

The agreement is excellent among the fits to three sample, and when the offsets are applied 
we obtain D*^ widths listed in the last row in Table The indicated uncertainties are only 
statistical. 
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Figure 1: Distribution of erg, the uncertainty on Q 
as determined from propagating track fitting errors. 
The arrow indicates a selection cut. 



Figure 2: Fit to nominal data sample. The different 
contributions to the fit are shown by different colors. 
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Figure 3: Fit to tracking selected data sample. 



Figure 4: Fit to kinematic selected data sample. 



We discuss the sources of systematic uncertainties on our measurements of the width of the 
D*^ in the order of their size. The most important contribution is the variation of the result as 
a function of the kinematic parameters of the D*^ decay. The next most important contribution 
comes from any mismodehng of ctq's dependence on the kinematic parameters. We take into 
account correlations among the less well measured parameters of the fit, such as fmis, and 
Cmis) by fixing each parameter at ±1 standard deviation from their central fit values, repeating 
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Figure 5: Likelihood function versus measured D*'^ width for the nominal (left), tracking (center), and kinematic 

(right) selected data samples. 



Table 3: Systematic uncertainties on the width of the D*+ and Qo 







Uncertainties in keV 










Sample 








Nominal 


Tracking 


Kinematic 


Source 


6r{D*+) 




6r{D*+) 








Running of Q 


16 


15 


16 


15 


16 


15 


Mismodeling of erg 


11 


< 1 


9 


4 


7 


< 1 


Fit Correlations 


8 


3 


9 


4 


9 


5 


Vertex Reconstruction 


4 


2 


4 


2 


4 


2 


Background Shape 


4 


< 1 


2 


< 1 


2 


< 1 


Offset Correction 


2 


NA 


6 


NA 


3 


NA 


Data Digitization 


1 


1 


1 


1 


1 


1 


Quadratic Sum 


22 


15 


22 


16 


20 


16 



the fit, and adding in quadrature the variation in the width of the D*'^ and Qq from their 
central values. We have studied in the simulation the sources of mismeasurement that give rise 
to the resolution on the width of the D*~^ by replacing the measured values with the generated 
values for various kinematic parameters of the decay products. We have then compared these 
uncertainties with analytic expressions for the uncertainties. The only source of resolution that 
we cannot account for in this way is a small distortion of the kinematics of the event caused by 
the algorithm used to reconstruct the origin point described above. We consider uncertainties 
from the background shape by allowing the coefficients of the background polynomial to float. 
Minor sources of uncertainty are from the width offsets derived from our simulation, and our 
data storage digitization resolution of 1 keV. 

An extra and dominant source of uncertainty on Qo is the energy scale of our measurements. 
We are still evaluating the size of this contribution. 

Table |3| summarizes the systematic uncertainties on the width of the D*^ and Qo- 
In summary we have measured the width of the D*^ by studying the distribution of the 
energy release in D*^ D^tt^ followed by — > K^i:'^ decay. With our estimate of the 
systematic uncertainties for each of the three samples being essentially the same we chose to 
report the result for the sample with the smallest statistical uncertainty, the minimally selected 
sample, and obtain 

r{D*+) = 96 ± 4 ± 22 keV, (3) 



where the first uncertainty is statistical and the second is systematic. 

This prehminary measurement is the first of the width of the D*"*", and corresponds to a 
strong coupling constants H 

5 = 0.59 ±0.01 ±0.07 or ^d-Dtt = 17.9 ± 0.3 ± 1.9. (4) 

This is consistent with theoretical predictions based on HQET and relativistic quark models, 
but higher than predictions based on QCD sum rules and lattice calculations. 
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